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SmBe was recently proposed to be both a strong topological insulator and a topological crystalline 
insulator. For this and related cubic topological Kondo insulators, we prove the existence of four 
different topological phases, distinguished by the sign of mirror Chern numbers. We characterize 
these phases in terms of simple observables, and we provide concrete tight-binding models for each 
phase. Based on theoretical and experimental results for SmBe we conclude that it realizes the phase 
with C(l" = -1-2, = -1-1, = —1, and we propose a corresponding minimal model. 


Topological insulators (TIs) with strong electronic cor¬ 
relations are considered to be of crucial importance in the 
exciting field of topological phases: They may provide TI 
states which are truly bulk-insulating - a property miss¬ 
ing from many Bi-based TIs - and they may host novel 
and yet unexplored interaction-driven phenomena. 

In this context, the material SmBg has attracted enor¬ 
mous attention: it has been proposed [1-3] to realize 
a three-dimensional (3D) topological Kondo insulator 
(TKI). This is a system where /-electron local moments 
form at intermediate temperatures T and are subse¬ 
quently screened at low T, such that a topologically non¬ 
trivial bandstructure emerges from Kondo screening [4]. 

While the results of numerous experiments on SmBe, 
such as transport [5-7] and quantum oscillation mea¬ 
surements [8] as well as angle-resolved photoemission 
spectroscopy (ARPES) [9-13] and spin-resolved ARPES 
(SP-ARPES) [14], appear consistent with the presence of 
Dirac-like surface states expected in a TKI, doubts have 
been raised about the proper interpretation of data [15- 
17]. Clearly, both experimental and theoretical progress 
is required for a thorough comprehension of this material. 

A recent theoretical insight [18] is that SmBg is also a 
topological crystalline insulator (TCI) [19], having three 
non-zero mirror Chern numbers. While Ref. 18 concluded 
that the latter are, modulo 2 or 4, independent of model 
details, we show in this Letter that the individual val¬ 
ues of these Chern numbers do depend on model details. 
As a result, we obtain four distinct TCI phases which in 
particular differ in their surface-state spin structure. Us¬ 
ing results from density-functional theory (DFT) and SP- 
ARPES we are able to single out one of the four phases 
as being relevant for SmBg. We also show that some of 
the previously employed models [3, 20] do not belong to 
this phase, and we propose a new minimal model for the 
correct TCI phase. Our results highlight the non-trivial 
role played by mirror Chern numbers in describing topo¬ 
logical band structures of actual materials. 

Lattice and symmetries. SmBg crystallizes in a simple- 
cubic (sc) structure, with its first Brillouin zone (BZ) 
shown in Fig. 1(a). In the following we use the lattice 
spacing a = 4.13 A as unit length. We may introduce mir¬ 


ror operators Mi = PC2{1), where P is the inversion and 
(72(0 is a rotation by tt about the axis I perpendicular to 
the mirror plane. For spin-1/2 particles Mf = — 1, such 
that Ml operators have eigenvalues ± i. For cubic symme¬ 
try there are three independent momentum-space planes 
which are invariant under the relevant mirror operators: 
the planes = 0 (equivalent to kx,y = 0) and fcz = tt 
(equivalent to kx,y = tt) are invariant under M^ while 
the plane k^ = ky (equivalent to k^ = —ky, ky = 
kz = ikx) is invariant under M^-y] see Fig. 1(a). 

Topological invariants. For each momentum-space 
mirror plane, one can define a mirror Chern number [18]: 
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with M|M^(k)) = ±i|u^(k)) and k lying in the plane 
BZ which is invariant under the symmetry operator M 
{M=Mz when BZ is = 0 or = tt, M = M^-y when 
BZ is kx = ky), and we sum over all N occupied bands. 
We note that -I- = 0 and, by cubic symmetry. 



Figure 1. (a) 3D BZ and its mirror planes kz = 0, 

kz = 'K, kx = ky-, high-symmetry points are F = (0,0,0), 
X = (7r, 0,0),(0,7r, 0), (0,0,7r), M = (0,7r,7r), (7r,0,7r), (7r,7r,0), 
R = (tTjTt, tt). (b) Corresponding 2D BZ for a (001) sur¬ 
face with its high-symmetry points F = (0,0), X = (7r,0), 
X' = (0, tt) and M = (tt, tt) and its mirror planes kx,y = 0, 
kx,y = TT, kx = diky, which are pairwise equivalent by Cav 
symmetry. (Planes at kx,y = tt and kx,y = —n are identi¬ 
cal and share a single label.) Also shown is the 2D Fermi 
surface of SmBe, consisting of three pockets centered at F, 
X, and X'. Dots mark its intersections with mirror planes. 
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Table I. Mirror Chern numbers ^q, _f, in different tight-binding models for SmBg, with INN (2NN) referring to 

first (second) neighbor hybridization. Also quoted are the number of bands Ni,ands, the orbital basis, and the non-zero terms 
in the ffamiltonian required by each model (divided into on-site energies e, d hopping f hopping , hybridization rf) [23]. 
The four phases, distinguished by u = sgn(Cj( ) and w = sgn(C(('^^QCj( are illustrated in Fig. 2 (where for Ft 

models* the spin directions needs to be reversed). 


= C'^ etc. A given insulating band- 

structure is thus characterized by a triplet of numbers 

(C£=m 

We recall that, according to bandstructure calculations 
[1, 2], SmBg is a strong TI: Band inversion between d 
and / states is achieved at the three X points of the 
3D BZ, leading to the parity invariants [24] (5(r) = -fl, 
(5(Ar) = —1, 6{M) = -1-1, S{R) = -fl and to the Z 2 in¬ 
dices {vq,viV2Vz) = (1,111). Ref. 18 proved, using the 
properties of the wavefunctions at high-symmetry points, 
that the mirror Chern numbers of SmBe obey =0 “ ^ 
mod 4, = 1 mod 4, = 1 mod 2 indepen¬ 

dent of microscopic details (like the Z 2 indices). Below we 
show this does not apply to the individual values of these 
Chern numbers which instead are model-dependent. 

Tight-binding models. Ab-initio bandstructure calcu¬ 
lations for SmBg [ 2 , 21 , 25] show that the bands close 
to the Fermi energy originate from Sm orbitals, namely 
in the d shell the Eg quartet \dx'^-y 2 f), \dx 2 -y 2 ),), 
1^22 f), 1^22 1 ) and in the / shell, where strong spin- 
orbit coupling pushes the j = 7/2 multiplet well above 
the Fermi energy, the j = 5/2 multiplet, which is 
split by cubic crystal field into a Fg quartet, [Fg^^i) = 

^1 ± I) + ^1 T I), |rf ±) = I ± i), and a Ft dou¬ 
blet, |r 7 ±) = -y/ll ± |) — T §)• Any reasonable 
model must thus consist of a sc lattice of Sm atoms, 
but the choice of included orbitals and tight-binding pa¬ 
rameters (on-site energies e and hopping parameters rj) is 
not univocal. In fact, published papers employed models 
with either 10 bands {Eg, Fg, Fy) [22], or 8 bands {Eg, 
Fg) [3, 22], or 6 bands {Eg, Fy) [22], or 4 bands (s dou- 
blet, Fg ' doublet) [20]. Once the basis orbitals are fixed, 
different sets of hopping parameters 77 may be used: in 
Table I we summarize some of these choices, with details 
given in the supplement [26]. 

We note that our discussion will be exclusively based 


on single-particle models: Although the Hubbard repul¬ 
sion among / electrons is not small, its effect at low tem¬ 
peratures can be captured in terms of renormalized /- 
electron kinetic energy and hybridization [27]. Hence, 
we can think of working directly with renormalized pa¬ 
rameters; our qualitative conclusions are independent of 
interaction-induced renormalizations. 

Distinct topological phases. Our main finding is that, 
upon computing mirror Chern numbers for different 
SmBg models, can be either -|-2 or — 2 and = 

-1-1 or —1, while = -fl always. This yields a total 

of four distinct phases, summarized in Table I. In particu¬ 
lar, models from Ref. 3 (here denoted as and Ref. 20 
(s/) are found to belong to the (—2, -|-1, —1) phase, while 
models from Ref. 22 with Fy and/or Fg multiplets real¬ 
ize either the (-1-2, -|-1, —1) or the (—2, -|-1, -|-1) phase; the 
(-1-2,-1-1,-1-1) phase is finally achieved in model ob¬ 
tained by modifying the hybridization term in model A^^ 
from Ref. 3. 

Phases and symmetries of surface states. In order to 
give a transparent physical meaning to the four phases, 
we relate the mirror Chern numbers to the properties of 
topological surface states. First, the absolute value of a 
mirror Chern number reflects the (minimum) number of 
Dirac points arising along a high-symmetry line in the 2D 
surface BZ [18, 20, 24]. For instance, on the (001) surface 
there are two Dirac points along the kx,y = 0 directions, 
and one along kx,y = tt and kx = Eky; see Fig. 1(b). 
This is consistent with the three Dirac cones at F, X 
and X' predicted by parity invariants [1-3] and observed 
experimentally [9-13]. Apparently, the absolute values of 
the C yield no additional information for the (001) surface 
(but they do for the (Oil) one [18]). This also shows that 
no other TCI phases are possible, since they would give 
rise to more Dirac cones than actually observed. 

However, the sign of a mirror Chern number gives new 
information, as it hxes the mirror-symmetry eigenvalue 
of surface states. For example, = -\-2 
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tells that along = 0 there are two bands with positive 
velocity along the zxx = y direction and eigenvalues 
+ i, while for = —2 these states have eigenvalue 

— i [20, 24], By a repeated use of this property we can 
assign all the ± i eigenvalues in Fig. 2(a)-(d) for the four 
possible phases. 

Surface-state spin structure. Since mirror eigenvalues 
are not directly measurable, we now link them to the spin 
expectation value (SEV) of surface states as observable 
in SP-ARPES experiments [14]. The idea is as follows: 
If each Dirac-cone state would simply behave as a spin- 
1/2 state with Mi = C 2 { 1 ) = = —iai, we would 

have Mi\si) = -i(Ti|sj) = -isi|s/) {I = x,y,x±y), with 
|s/) an eigenvector of the Pauli matrix ai with eigenvalue 
Si = ±1. Hence, Mi eigenvalues of ± i would immediately 
give the direction of the spin, e.g., an M^ eigenvalue of 

— i would imply a spin pointing towards -\-x. However, 

additional minus signs arise for orbitals odd under M;, 
and, due to the spin-orbit coupling in the / shell, surface 
states are not eigenstates of spin operators. By explicitly 
invoking the properties of the basis states we can de¬ 
rive rules relating the mirror eigenvalues to the SEV: We 
find that M;|'0) = — i ['i/') implies that the SEV for state 
I'l/') points along direction ±/, with the plus sign always 
holding except for Mx±y acting on and Pg^^ states 

and for M^^y on Py states [26]. Using these rules, and the 
fact that for small momenta around P the SEV does not 
depend on the eigenvalues of Mx±y (see below), we can 
construct a full portrait of the SEV for each phase once 
a basis has been fixed, see Fig. 2 for Pg states. 

We proceed by analyzing these portraits in detail. In¬ 
terpreting the SEV in Fig. 2 as true spin [28], the effective 
low-energy model at the X points can be written as: 

Hx = vi{kx - TT)ay - V2kyaoc, ( 2 ) 

Hx'= V2kxcry - vi{ky - TT)ax. ( 3 ) 

Importantly, the velocities vi ^2 have the same (opposite) 
sign if C'^ > 0 (< 0), respectively. We introduce 

w = sgn(C/'^^) = sgn(uiU 2 ) which describes how 
the mirror eigenvalues and the SEV on the X cones evolve 
upon 7r/2 momentum-space rotations around X. This 
implies that w is a winding number, when rc = -|-1 spin 
and momentum rotate around X in the same direction. 
Figs. 2(b,d), while for w = — 1 they rotate in the opposite 
direction. Figs. 2(a,c). We note that w = —1 is only 
allowed for cones at low-symmetry points, see below. 

Similarly, we introduce v = sgn(C'//^gC//^^^) which 
now describes how the mirror eigenvalues on the F cone 
evolve upon 7r/4 rotations. It turns out, however, that 
V does not affect the SEV of the f cone: For sufficiently 
small momenta the system must display continuous ro¬ 
tation symmetry around F which implies a fixed angle 
between SEV and momentum along the contour (i.e., 
a winding number -pl). Thus the effective low-energy 



Figure 2. Surface-state mirror eigenvalues ±i (defined at 
snrface momenta where the Fermi contour crosses a mir¬ 
ror plane) for the four phases with mirror Chern numbers 
(^fc) =0 > =7r I ) as quoted in the panels. The C signs 

can be condensed into v = sgn(C/'^^gCj!'^^j,^) and w = 
sgn(C)/^.„.), see text. Green arrows show the surface- 
state SEV for Fg states; for r7 states (see Table I for models) 
it would be reversed. In the FTB model this reversal depends 
on the relative weight of Fy and Fg states, with Fg dominating 
for realistic parameters [22]. 

model at the F point is independent of u, and reads 

iFp = vow{kx(Ty - kyax) (vo > 0), (4) 

with spin and momentum always forming a mutual angle 
±7r/2 given by the sign of vqw, i.e., the chirality. We 
note that the w factor in Up is needed to give the cor¬ 
rect relative spin direction along k^ and ky between the 
f and the X cones, which is fixed by jC/’^gl = 2. Using 
Up we can show [26] that v instead dictates symmetry 
properties of the states composing the F cone for small 
momenta: for u = -pl these are (i ^2 and Fg ^ (Xg symme¬ 
try representation) whereas for v = —1 these are 
Fg^^ and Fy (corresponding to Xy). 

Relation between v,w and electronic structure. In order 
to determine the TCI phase of SmBe we now connect the 
parameters v and w to the bulk bandstructure. Surface 
states near F can be computed by perturbatively expand¬ 
ing the Hamiltonian around X = (0,0, tt) [29, 30], where 
it decouples into a subspace with symmetry representa¬ 
tion Xy, composed from dx^-y^, Fg^^ and Fy states, and 
a subspace with symmetry representation Xg and d^^, 
Fg ' states. Band inversion can be achieved in only one 
of these subspaces: in model from Refs. 3 and 30 it 
is assumed to be in the subspace Xg, while DFT calcu- 
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lations [25] display band inversion in subspace X^j. Since 
surface states at f only exist in the subspace where band 
inversion is achieved [30], and this depends on which d 
band has a minimum at X, we can link these two options 
to the symmetry of the d-band minimum at X: {Xq ) 

in the first case, leading to v = + 1 , dj; 2 -y 2 (X^) in the 
second, leading to u = — 1. Since we expect DFT to be 
reliable for weakly correlated orbitals, we conclude that 
u = — 1 in real SmBg. 

Once the band-inversion subspace, and thus v, is fixed 
by the choice of the hopping terms, w depends on the 
hybridization term: out of those which lead to a full gap, 
some lead to w = - 1 - 1 , some others to rc = — 1 , as shown 
in Table I. When considering Tg states, DFT calculations 
[ 21 , 22 ] show the largest of the hybridization terms to 
be [23], which, however, alone does not lead to a 
gap for u = — 1 , since it does not couple the inverted 
bands d^^-y^ and Fg^^ along F-X for symmetry reasons. 
The second-largest term is which gives re = -fl (r|^ 
in Table I). However, a competition with the tc = — 1 
phase can be observed when retaining both Fy and Fg 
multiplets [ 22 ], since the Fy doublet alone with Eg is in 
the w = —1 phase, ryy^ being the largest hybridization 
term in this case (Fy^). 

Choice of phase and minimal model. We now make 
use of the following ingredients: (i) The band inversion 
properties from DFT calculations yield v = —1. (ii) The 
SP-ARPES experiment of Ref. 14 shows a winding num¬ 
ber on X cones oiw = -1-1. Together, this uniquely yields 
(-1-2, -1-1, —1) as the best candidate TCI phase for SmBg. 

At the same time, we realize that the model, ob¬ 
tained retaining the Eg and Fg quartets, with ~ 

~ —3/8 and the 77 /^ hybridization term, is 
the simplest one giving the correct mirror Chern num¬ 
bers (- 1 - 2 ,-|- 1 , — 1 ) and reproducing as closely as possi¬ 
ble the DFT bandstructure. We propose it as a new 
minimal model for SmBg, with a concrete tight-binding 
parametrization given in the supplement [26]. We stress 
that this conclusion is based on the experimental result 
w = -|-l;w = —1 would instead lead to a minimal model 
entailing the Eg quartet and the Fy doublet (Fy^). The 
competition between the w = - 1-1 and w = —1 phases, 
emerging from the Fg quartet and the Fy doublet, respec¬ 
tively, depends on model details which have not been ob¬ 
tained with sufficient accuracy from DFT. This applies 
for instance to the Fy-Fg energy difference whose vari¬ 
ation can drive a topological phase transition between 
states with different w. 

Quasiparticle interference (QPI). QPI can be used as 
a probe of the topological character of surface states [31- 
33], owing to their spin texture; in TCIs, impurities which 
do not break mirror symmetries cannot induce transi¬ 
tions between states with opposite mirror eigenvalues 
[34]. Here we show how the winding number w affects 
intercone scattering. Comparing the spin portraits for 
two relevant models in Figs. 3(a) and (c) we see that, for 



Figure 3. Isoenergy contours with SEV (a,c) and QPI signal 
(b,d) at the Fermi level for two different SmBe tight-binding 
models. Large arrows show intercone scattering wavevectors. 
(a,b): Four-band model sf in the (—2,-|-l,—1) phase (w = 
— 1); (c,d): Eight-band model Fg^ in the (-1-2,-|-1, —1) phase 
(w = -t-1); with parameters given in the supplement [26]. An 
intercone QPI signal is seen in (b), but not in (d), dictated 
by the surface-state spin structure. 


pairs of stationary points [35] (states which have the same 
tangent to the isoenergy contour) belonging to different 
X cones, w = —1 (w = -1-1) implies roughly parallel 
(antiparallel) spin, respectively. Consequently, w = —1 
yields a QPI peak from intercone scattering [36, 37], 
Fig. 3(b), while w = -1-1 does not [22], Fig. 3(d). Thus, 
QPI can probe w by looking at intercone scattering, and 
hence provides information on the spin structure of the 
Dirac cones and, indirectly, on mirror Chern numbers. 

Conclusions. In this Letter we have studied mirror 
Chern numbers in different models for the correlated 
topological insulator SmBg. The absolute values of the 
mirror Chern numbers are fixed by symmetry proper¬ 
ties of the atomic orbitals and the bandstructure and 
are thus independent of model details, while their signs 
are model-dependent. This yields four possible phases 
- all sharing the Z 2 topological indices (1,111) - which 
can be characterized by the two combinations of signs 
V = sgn(C+^^oC+^fcJ and w = sgn{C+^^Cl^J. These 
have a straightforward physical interpretation in terms 
of the symmetry of the d-band minimum at X and 
the spin texture of the X cone on a (001) surface, re¬ 
spectively. Our analysis constrains SmBg to be in the 
iCt=o^^t=^^^t=kJ = (+2,+1,-1) phase, by using 
the results of DFT calculations and of SP-ARPES [14] 
which give v = —1 and rc = -1-1, respectively. Eor the 
Tu = -|-1 case we predict a weak QPI signal from inter¬ 
cone scattering. 

From our analysis we propose a minimal tight-binding 
model for SmBg consisting of 8 orbitals (the Eg and Fg 
quartets). The model equals the one of Ref. 3 in the 
choice of the orbital basis but differs in both the kinetic 
energy and the hybridization parameters; it better re¬ 
produces DFT and SP-ARPES data, at the price of in¬ 
troducing second-neighbor hopping parameters. We fi¬ 
nally speculate that chemical substitution can tune both 
the Fr-Fg energy difference and the hopping parameters, 
possibly leading to a topological phase transition from 
the ic = -1-1 to the w = — 1 state. Such a transition is 
accompanied by the closing of the bulk gap. 
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We believe that our results advance the understanding 
of the complex TKI material SmBg (and its cousin PuBg 
[21]), and potentially of other TIs, too. 

This research was supported by the DFG through SFB 
1143 and GRK 1621 as well as by the Helmholtz associ¬ 
ation through VI-521. 

Note added: At the day of submission of this paper, 
a related paper [38] appeared on arXiv, which reaches 
conclusions compatible with ours. For comparison we 
note that their Ghern number Cd equals our {—C^ ). 
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I. MODEL HAMILTONIANS 

In this section we provide explicit expressions for the 
tight-binding Hamiltonians referred to in the main text 
and Table 1 . 


A. Parametrization 

We will use a, r, and ^ Pauli matrices as follows: a 
matrices act in the (pseudo) spin subspace, r in the or¬ 
bital subspace, and /x in the d-f subspace. Then, terms 
proportional to /Xq and /x^ denote the kinetic energy (with 


subscript 0 representing the 2 x 2 unit matrix) while 
terms proportional to /Xj, correspond to d-f hybridiza¬ 
tion. For convenience, we introduce /xi = (/xq -f /Xz)/ 2 , 
d 2 = (mo - Mz)/ 2 , ti = (to -I- rz)/ 2 , T 2 = (tq - t^)/ 2 : yxi 
projects onto the d subspace, /X2 onto the / subspace, ti 
onto the orbital subspace 1 {d,^2_y2, Pg^^), and T2 onto 
the orbital subspace 2 (d^z, Pg ). We also employ abbre¬ 
viations Cx = coskx, Cy = cosky, Cz = cos kz, Sx = sinkx, 
Sy = sin ky, Sz = sinfcz. 


We consider three classes of models: 




8 -band models consisting of the Eg{d) quartet interacting with the r 8 (/) quartet. 
The parameters are = {dx2^y2\H\dx2^y2)ooo = {dz2\H\dz2)ooo, = 
(Pf^ + |iJ|rf 4 )ooo, = {dx2.y2\H\dx2.y2)ooi, rjf = {dz2\H\dz2)ooi, 

r;/! = (PW + |id|r(i4)ooi, r;/i = (p(2) + |id|r(24)ooi, V 

Tlf = {dx2.y2 t |iP|r^' 4 )ooi, Vf = {d.2 t |dd|pf+)ooi, Vf 

= {dz2 t |id|P 8 ^^-f)iio, and Hamiltonian H = He^ + Hug + VEg-rs- 


(pw + |dd|r('4)ooo 

r,f = {dz2\H\dz2) 

+ |id|rf+) 

= {dx2.y2 t \H\T‘i^+) 


110? 

110 ? 

110 ? 


Hifk) = e^'Toaotii - ttaoHi ■ 

■{[{Cx +Cy-\- Cz){rfz + v'x) + ‘^dfi.CxCy + CyCz + CzCx)]to -b 
+ [(cx + Cy- 2cz)(jtz - Vx) + 2??f (-2cg;Cy -b CyCz -b C^;)]T^/2 -b 

-b[vl5(c,, - Cy){T]l^ - + 2cz'nf)]Tx/2}, (SI) 

with i = Eg{d) = 1,P8(/) = 2, 

V£;,-rs(k) = -^iyV{sxCrxMrif -b rjf + 2{cz + Cy)(jqf + 77 ^)] -b 
+SyayTo[Tif -b ryf + 2 (c^ -b Cx){vf + vf)] + 

+SzCr;,To[? 7 "^ -b Tif + 2{cx + Cy){r)f + r]f)] + 

+sxaxrz[-v:^ + of + 2(c. - 2cy){rif - ryf )]/2 + 

+Sy(JyTz[-rif + of + fcz - 2cx){vf - )]/2 + 

+SzCrzTz[r]f - of + fx + Cy){r]f - r/f)] -b 

+fisxaxTx[of - of + ^cziof - df )]/2 + 

+V3syayTx[-of + of - ‘^Cziof - )]/2} + 

+ViszazTx{cx-Cy){r]f-ijf)}; (S2) 

• 6 -band models consisting of the Eg{d) quartet interacting with the r 7 (/) doublet. 

Here the parameters are e"*, e( = (P 7 -b |H|r 7 +)ooo, of, of, of, of = (^7 + \H\'^ 7 +)ooi, 
Tjf = (P 7 + |id|r7+)iio, of = (r7 + |H|r7+)iii, of = {dx^-y^ t |i^|r7+)ooi, of = (4= t |i?|r7+)iio; 
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and Hamiltonian H = He^ + 

Hrjik) = -t^^(7o^i2[2v7^{cx + Cy + c^) + 

+4??7^(Ca;Cy + Cj/C^ + CzCx) + , (S3) 

VEg_r7(k) = -fj.yV{2s^a:,uJo[{r]^^ - V3(c^ + Cy)rif] + 

~^{sxO'x + SyO'y)uJo[—ri^^ + 2^/3czT]^‘^] + 

+2s2(JzWi(Cj, - Cx)?77^ + 

+Sa;(Ja;a;i [—+ 2r]^‘^(^Cz + 2cy)] + 

+Sy(jyWi[v^?77^ - 2rif{cz + 2cx)]}, (S4) 

where wq = (1,0) and wi = (0,1) act in the orbital space, projecting respectively onto the d^^-y^ and dz 2 
doublets, and He^ (k) is the same as from the 8-band model; 

• 4-band model consisting of conduction (s) electrons interacting with the rg^^(/) doublet. 

This has parameters = (s|i?|s)ooo, e-^ = -b |iJ|rf^-b)ooo, = (s|-ff|s)ooi, = (s|H|s)iio, 

= (rf ^ -b |i?|rf^-b)ooi, = (rf ^ + |H|rf^-b)iio, = (s t l-ff|r8^^-b)ooi: and Hamiltonian 

H = Hg + H (2) -b V „(2): 

^8 * 8 

Hi{k) = -2ti^i,[Tf^{cx + Cy + Cz) + 2r]"‘^{CxCy + CyCz + CzCx)]cro, (S5) 

V(,_p(2)(k) = 2V fiy-q'"^ {Sxdx + Syay + SzCTz), (S6) 

with i = s= l,r^g\f) = 2. 


We now relate these general model Hamiltonians to 
models which have appeared in the literature and to those 
which are listed in Table 1 of the main text. 

• The model by Legner et aid (here abbreviated 
as sf) is identical to the 4-band model, setting 
r^s2jrjsi „ rjf ^^ —0.5. 

• The 8 -band model by Alexandrov et aid (A^^) is 

obtained in the Eg-Ts basis by setting ~ 

~ —0.3, and keeping the 77 "^ INN term 
in the hybridization. For completeness, we have 
also considered the same model with the 77 “^ (A^“) 
2NN term in the hybridization, since it belongs to 
a different phase. 

• Our Tg^ and r|^ 8 -band models are obtained re¬ 
taining the Eg and the Tg quartets, with 77 ^^/ 772 ^ ~ 
—3/8, r][‘^/r]E ^ —1/2 and keeping either the INN 
77 "^ (Tg^) or the 2NN 77 /^ (r|^) term in the hy¬ 
bridization. 

• Our and T^^ 6 -band models are obtained re¬ 
taining the Eg quartet and the r 7 doublet, with 

vt'^/Vz^ ~ —3/8, 777 ^/ 777 ^ ~ 1/2, and keeping ei¬ 
ther the INN 777 ^ (r^^^) or the 2NN 777 ^ (Ff^) term 
in the hybridization. 


• The full tight-binding calculation [FTB) of Ref. 3 
comprises the Eg and the Fg quartets, and the r 7 
doublet, for a total of 10 orbitals per site. Param¬ 
eters are retained up to 7NN, such that we cannot 
give a complete analytic expression. 

• The model by Takimoto^ (T) is equivalent to our 
Eg — model, retaining, however, in addition to 
more / terms, also the 77 /^ hybridization term. This 
term alone, in our approach, does not lead to an 
insulator, since it does not provide a gap along the 
F — A direction in the BZ. Therefore we have not 
discussed this model in the paper. 

B. Realistic tight-binding parameters 

A plausible parameter set for the F|^ model, which 
essentially reproduces the DFT bandstructure of PuBg,^ 
is:^ td = leV, t/ = 0.01 eV, Cg — = —2eV, 77 /^ = 

0.8, 77 f = -0.3, 77/1 = -4, 77/2 = 2, P 77/2 = 0.06 eV. 
Similarly, a plausible parameter set for the F^^ model is: 
td = leSf, tf = 0.01 eV, e( - Sd = -1.9eV, r]f = 0.8, 
r,f = -0.3, r](^ = 2.5, 77 /^ = 1.25, ^ 77^2 o.05eV. 

The effect of electronic correlations, mainly absent from 
DFT, will be a downward renormalization of t/ and V, 
see Ref. 3 for a more extensive discussion. 
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Panels (a) and (b) of Fig. 3 of the main text have 
been calculated using the sf model with ts = leV, 
tf = O.OOSeV, ef - €2 = -2eV, = 1, 

j^s 2 _ ^/2 _ _Q 5 ^ Vr]^^ = 0.2 eV. Panels (c) and (d) of 
the same figure have been obtained using the r|^ model 
with td = I eV, tf =0.01 eV, Cg — = —2 eV, = 0.8, 

ryf = -0.3, ryf = -4, = 2, Vrjf = 0.02 eV, 

Vr]l^ = —O.OGeV. Compared to the parameter set from 
above, this includes the additional hybridization param¬ 
eter which improves the quantitative agreement with 
DFT, but does not change the topological properties. 

II. SEV FROM MIRROR SYMMETRY 
EIGENVALUES 

Here we show how to compute the spin expectation 
value (SEV) of surface states from the corresponding mir¬ 
ror eigenvalues. Mirror operators are defined as Mi = 
where I ^ -I, {I = x,y,z), : 

X O =Fy, and = —iai {I = x,y,z), 

^x±y = - ± cry)/V2 = - i a^±y. 

We now determine the action of these mirror operators 
on Eg, Fy, Fg states: 

Mi\dmcr) = m = x"^ - j/^, 2 :^; cr =t,i, (S7) 

Mi|r„a) = +ia,|r„(T), m = 8(i),8(2),7;a = ±, (S8) 


where I = x,y, z, and ct/ acts on the subspace spanned 
by t, i for d states, or by -|-, — for / states; moreover 


^x:k.y\^x'^ “1“^ \/2 \^x‘^—y‘^^') "> 

(S9) 

Mx±y\dz2a) = 1 ^ ^|dz2Cr), 

(SIO) 

Mx±y\Tm<y)= 1 ^ ^IF^ct), m 

= 8(^\7,(S11) 

A7x±y|r„(T) =- 1-1 ^ ^IF^ct), m 

= 8(2). (S12) 

We can now find which states belong to eigenvalues 
(±i). We start with M^- States with eigenvalue (—i) 

are: 


IV’Dx) = Mx2-y2)(l, 1)/V2, 

(S13) 

IV’2■:x) = M.= )(l;l)/^/2, 

(S14) 

iv’3:x) = i4'^)(1;-i)/V2, 

(S15) 

IV’4:x) = |rr)(i,-i)/V2, 

(S16) 

IV’5■:x) = l^7)(l,-l)/^/2 

(S17) 


where the two-dimensional vector is in the (f, i) {d 
states) or (-I-, —) (/ states) subspace. The states IV’fa;) 
and \' 4 ’ 2 x) are spin eigenstates, pointing towards the 
positive X (-fx) direction, but the other three states 
are not, because -I- and — are not spin indices. We 
can, however, compute the SEV on these states with 


the use of Table SI. Now, on Fg states, (ctj,, (Jy, ct^) ~ 
{—Sx, ~Sy, —Sz), where on the left-hand side at are true 
spin operators, while on the right-hand side Si act on the 
pseudospin space; for Fy states, instead, {o'x,o'y,az) ^ 
(+ 5 x 1 +Sy, TSz). So, for states and IV '4 x) the SEV 

is pointing in the +x direction, while for IV'^^,) it is point¬ 
ing in the —x direction. 

To conclude, for a state with Mx|V’) = — i |^)j even 
though it is not in general an eigenstate of the spin, its 
SEV is pointing along -\-x if it is a d or a Fg state, but 
along —X if it is a Fy state. The opposite holds for eigen¬ 
value -I- i. We can repeat the same reasoning for My, and 
conclude that, for a state \\p) with My\xp) = — i|V') hs 
SEV is pointing along +y if it is a d or a Fg state, but 
along —y if it is a Fy state. 

For Mx-y the eigenstates with eigenvalue (— i) are: 

IV'px-y) = (SI 8 ) 

IV' 2 ■:x-,) = l4^)(l+■'"/")/^/2, (S19) 

IV'3■:x-,) = l4'^)(l+-‘"/")/^/2, (S20) 

IV'W = (S21) 

IV'5':x-,) = |r7)(l,e-‘"/")/y2, (S22) 

which means that |V'i)x-y) ^ ®Pi^ eigenstate with spin 

pointing along (—x -I- y), while {ip^x-y) points along 
(x — y); for IV’g^x-y) the SEV points along (—x -I- y), for 
\'^ 4 ,x-y) along (x - y), and for along (x - y). 

From this we read off that for a state |^/>) with 
Mx_y|f/') = — i|V’) the SEV points along (x — y) if it 
is dz 2 , Fg ' or Fy, and along (—x -I- y) if it is dx^-y^ or 
Fg^l. Similarly, for a state |'(/>) with Mx+ylV’) = ~i|V’) 
the SEV points along (x -I- yy) if it is d 22 , Fg or Fy, and 
along (—X — y) if it is dx2-y2 or Fg^l. This concludes the 
set of rules used in the main text. 

( 2 ) 

In addition, the 4-band model requires s and Fg 
states: s states behave as dz 2 states, so if Mx\ip) = 
— ilV'), the spin points along +x, if Mylijj) = —i !■!/:), 
it points along +y, and if Mx-y\ip) = — i|f/'), the spin 

points along (-|-x —jy). As for Fg , if Mx\ip) = — i lip), the 
SEV points along +x, if My\\p) = — i | 7 />), it points along 
+y, and if Mx_y IV’) = — i IV’); + points along (-|-x — y). 
We note that the transformation x —> —x, y —)■ —y, 
needed to obtain the hybridization term Eq. (S 6 ) with 
all plus signs, changes the explicit form of Fg states, 
but neither their symmetries properties nor their SEV. 


III. SYMMETRY OF SURFACE STATES AT P 

In this section we show how v determines the symmetry 
of the states which build the F cone for small momenta. 

As usual® we introduce a time-reversal doublet IV’-i-); 
IV’-) which describes surfaces states at f; it can be writ- 
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Table SI. SEV in the Ty - Tg basis. 


ten as: 


This state has the properties 


IV'S-) — t) + 4>2\dz2 f) + 

T'/'sirg ^+) + </!'4|rg ^+) + ^5|r7+), (S23) 
IV’-) = <t>l\dz: 2 _y 2 1 ) + (I) 2 \dz 2 i) - 

-Valr^-) - -) - Vslrr-). (S24) 

(S25) 

It is constructed such that 


M,|0+(O))=-i|V.+ (O)), 

A4‘a;|V)+(7r/2)) =+i|(/)+(7r/2)), 
Mz:-y\(j)''~{TT/A)) = - i |V)+(7r/4)), 
M^_y|(/)+(7r/4)) =+i|(/)+(7r/4)), 


(529) 

(530) 

02, <^4 = 0, (S31) 

01,03,05 =0.(S32) 


((cTa;, (Jy, az)} = {aSx, asy, bsz), (S26) 

with a and b real numbers, as can be seen through the use 
of Table SI, and where we suppose that the Tg coefficients 
03, 04 and the Fy one 05 are not non-zero at the same 
time. The eigenstate with positive energy for the Dirac 
Hamiltonian iFp = vow{kxO'y — kyU^) with ru = 4-1 is 
now 

l</'+(^^k)) = ^(|0+)+ie'^-‘|0_)), (S27) 

cosdk = ^, sindk = ^- (S28) 

k k 


We now see that when ■u; = -l-l,u = -|-l, we must have 
Ma;_j^|0+(7r/4)) = -bi 10+( 75 / 4 )), so, in this case, the 

r cone is entirely composed by dz 2 and Tg ; however, 
when w = 4-1, v = —1, we must have |0+(7r/4)) = 
— i 10+( 77 / 4 )), so the f cone is now entirely composed by 
dx 2 -y 2 , Tg^^ and Fy. For w = — 1 all mirror eigenvalue 
signs and spin directions are reversed, but the results are 
the same. We conclude that v dictates the symmetry of 
the surface states at f for small momenta. 
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